Ultra-High Energy Neutrino Fluxes and Their Constraints 
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Applying our recently developed propagation code we review extragalactic neutrino fluxes above 
10 14 eV in various scenarios and how they are constrained by current data. We specifically identify 
scenarios in which the cosmogenic neutrino flux above ~ 10 18 eV, produced by pion production 
of ultra high energy cosmic rays outside their sources, is considerably higher than the "Waxman- 
Bahcall bound" . This is easy to achieve for sources with hard injection spectra and luminosities that 
were higher in the past. Such fluxes would significantly increase the chances to detect ultra- high 
energy neutrinos with experiments currently under construction or in the proposal stage. 



I. INTRODUCTION 



The farthest source so far observed in neutrinos was su- 
pernova SN1987A from which about 20 neutrinos in the 
10-40 MeV range where detected JD. Extraterrestrial 
neutrinos of much higher energy are usually expected 
to be produced as secondaries of cosmic rays interact- 
ing with ambient matter and photon fields and should 
thus be associated with cosmic ray sources ranging from 
our Galaxy to powerful active galactic nuclei (AGN) ||. 
Traditional neutrino telescopes now under construction 
aim to detect such neutrinos up to ~ 10 16 eV by look- 
ing for showers and/or tracks from charged leptons pro- 
duced by charged current reactions of neutrinos in ice, 
in case of AMANDA j| |J] and its next generation ver- 
sion ICECUBE [§, in water, in case of BAIKAL §, |, 
ANTARES §, and NESTOR 0, or underground, in 
case of MACRO ]ljj (for recent reviews of neutrino tele- 
scopes see Ref. [|12|). 

On the other hand, the problem of the origin of the cos- 
mic rays themselves is still unsolved, especially at ultra- 
high energies (UHE) above ~ 4 x 10 19 eV, where they 
lose energy rapidly by pion production and pair pro- 
duction (protons only) on the cosmic microwave back- 

For sources further away than 
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a few dozen Mpc this would predict a break in the cos- 
mic ray flux known as Greisen-Zatsepin-Kuzmin (GZK) 
cutoff pq |, around 50EeV. This break has not been 
observed by experiments such as Fly's Eye [jl6| , Hav- 
crah Park 0, Yakutsk fl§ and AGASA @, which in- 
stead show an extension beyond the expected GZK cutoff 
and events above 100 EeV (however, the new experiment 
HiRes ]2Cj currently seems to see a cutoff in the monoc- 
ular data ^l)). This has lead to the current construc- 
tion of the southern site of the Pierre Auger Observa- 
tory [p2[ , a combination of an array of charged particle 
detectors with fluorescence telescopes for air showers pro- 
duced by cosmic rays above ~ 10 19 eV which will lead to 
an about hundred-fold increase of data. The telescope 
array, another planned project based on the fluorescence 



technique, may serve as the optical component of the 
northern Pierre Auger site planned for the future p3| . 
There are also plans for space based observatories such 
as EUSO Q and OWL [^7| of even bigger acceptance. 
These instruments will also have considerable sensitivity 
to neutrinos above ~ 10 19 eV, typically from the near- 
horizontal air-showers that are produced by them |p8| . 
Furthermore, the old Fly's Eye experiment |2^] and the 
AGASA experiment |3(| have established upper limits on 
neutrino fluxes based on the non-observation of horizon- 
tal air showers. 

In addition, there are plans to construct telescopes 
to detect fluorescence and Cerenkov light from near- 
horizontal showers produced in mountain targets by neu- 
trinos in the intermediate window of energies between 
- 10 15 eV and - 10 19 eV || H- The alternative of 
detecting neutrinos by triggering onto the radio pulses 
from neutrino-induced air showers is also investigated 
currently |35| ]. Two implementations of this technique, 
RICE, a small array of radio antennas in the South pole 
ice |34j , and the Goldstone Lunar Ultra-high energy neu- 
trino Experiment (GLUE), based on monitoring of the 
moons rim with the NASA Goldstone radio telescope for 
radio pulses from neutrino- induced showers p5[ , have so 
far produced neutrino flux upper limits. Acoustic detec- 
tion of neutrino induced interactions is also being con- 
sidered 

The neutrino detection rates for all future instruments 
will crucially depend on the fluxes expected in various 
scenarios. The flux of "cosmogenic" neutrinos created 
by primary protons above the GZK cutoff in interac- 
tions with CMB photons depends both on primary pro- 
ton spectrum and on the location of the sources. The 
cosmogenic neutrino flux is the only one that is guaran- 
teed to exist just by the observations of ultra-high energy 
cosmic rays (UHECRs) and was studied soon after the 
discovery of the CMB |j37j . Note, however, that there 
is no firm lower bound on the cosmogenic neutrino flux 
if the UHECR sources are much closer than the GZK 
distance. 

If sources are located beyond the GZK distance and the 
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proton flux extends beyond the GZK cutoff, the neutrino 
fluxes can be significant. This possibility is favored by the 
lack of nearby sources and by the hardening of the cosmic 
ray spectrum above the "ankle" at ~ 5 x 1 18 eV. It is 
also suggested by possible correlations of UHECR arrival 
directions with compact radio quasars [j38| and more sig- 
nificant correlations with BL Lacertae objects [p9|, some 
of which possibly also luminous in GeV 7— rays pCfl and 
at distances too large to be consistent with the absence 
of the GZK cutoff. 

Whereas roughly homogeneously distributed proton 
sources can naturally explain the UHECR flux below the 
GZK cutoff (see, e.g., Refs. j|l| and g|), the highest en- 
ergy events may represent a new component. They may 
be new messenger particles, which propagate through the 
Universe without interacting with the CMB E3] , or may 
have originated as extremely high energy {E~£ 10 23 eV) 
photons, which can propagate several hundred Mpc (con- 
stantly loosing energy) and can create secondary photons 
inside the GZK volume [[l2| . Decaying super-heavy relics 
from the early Universe (see Ref. H, H for reviews) can 
also explain UHECRs and predict UHE neutrino fluxes 
detectable by future experiments. 

Another speculative possibility is to explain UHECRs 
beyond the GZK cutoff by the UHE protons and photons 
from decaying Z-bosons produced by UHE neutrinos in- 
teracting with the relic neutrino background The 
big drawback of this scenario is the need of enormous pri- 
mary neutrino fluxes that cannot be produced by known 
astrophysical acceleration sources without overproducing 
the GeV photon background ^6|, and thus most likely 
requires a more exotic top-down type source such as X 
particles exclusively decaying into neutrinos |47| . As will 
be shown in Sect. VI, even this possibility is also signifi- 
cantly constrained by existing measurements. 

Active galactic nuclei (AGN) can be sources of neutri- 
nos if protons are accelerated in them . In the present 
paper we consider only the two representative limits of 
low and high optical depth for pion (and neutrino) pro- 
duction in the source. In the first case the protons accel- 
erated in the AGN freely escape and neutrinos are pro- 
duced only in interactions with the CMB (cosmogenic 
neutrinos). For the second case we discuss an example of 
possible high neutrino fluxes from a non-shock accelera- 
tion AGN model [Q, in which primary protons lose all 
their energy and produce neutrinos directly in the AGN 
core. 

Motivated by the increased experimental prospects for 
ultra-high energy neutrino detection, in the present paper 
we reconsider flux predictions in the above scenarios with 
our recently combined propagation codes Jl2[ [f(], ^lj . 
Our main emphasize is thereby on model independent 
flux ranges consistent with all present data on cosmic 
and 7— rays. For any scenario involving pion production 
the fluences of the latter are comparable to the neutrino 
fluences. However, electromagnetic (EM) energy injected 
above ~ 10 15 eV cascades down to below the pair pro- 
duction threshold for photons on the CMB and EM en- 



ergy above 100 GeV also cascades down due to the pair 
production on infrared/optical background. The result- 
ing intensity and spectrum of 7— rays below 100 GeV is 
rather insensitive to these backgrounds |3^, HH (for re- 
view see 14 ). The cascade thus gives rise to a diffuse 
photon flux in the GeV range which is constrained by 
the flux observed by the EGRET instrument on board 
the Compton 7— ray observatory pif . For all neutrino 
flux scenarios the related 7— ray and cosmic ray fluxes 
have to be consistent with the EGRET and cosmic ray 
data, respectively. 

Sect. II summarizes the numerical technique used in 
this paper. In Sect. Ill we discuss the cosmogenic neu- 
trino flux and its dependence on various source charac- 
teristics. We specifically find an upper limit considerably 
higher than typical fluxes in the literature and remark 
why it is higher than the Waxman-Bahcall (WB) | 5q] 
and even the Mannheim-Protheroe-Rachen (MPR) |5(| 
bounds for sources transparent to cosmic and 7— rays. 
In Sect. IV we review neutrino flux predictions in top 
down scenarios where UHECRs are produced in decays of 
super-massive particles continuously released from topo- 
logical defect relics from the early Universe. Sect. V 
discusses neutrino fluxes in scenarios where the cosmic 
rays observed at the highest energies are produced as 
secondaries of these neutrinos from interactions with the 
relic cosmological neutrino background, often called Z- 
burst scenario. In Sect. VI we focus on a combina- 
tion of top-down and Z-burst scenarios considered by 
Gelmini and Kusenko p7| , namely super- heavy parti- 
cles mono-energetically decaying exclusively into neutri- 
nos" (see, however, (48)). In Sect. VII we discuss possible 
high neutrino fluxes from a non-shock acceleration AGN 
model 1 49 . Finally, in Sect. VIII we conclude. 



II. NUMERICAL TECHNIQUE 

Our simulations are based on two independent codes 
that have extensively been compared down to the level of 
individual interactions. Both of them are implicit trans- 
port codes that evolve the spectra of nucleons, 7— rays, 
electrons, electron-, muon-, and tau-neutrinos, and their 
antiparticles along straight lines. Arbitrary injection 
spectra and redshift distributions can be specified for 
the sources and all relevant strong, electromagnetic, and 
weak interactions have been implemented. For details see 
Refs. || 013. 

Relevant neutrino interactions for the Z-burst scenario 
are both the s-channel production of Z bosons and the 
t-channel production of W bosons. The decay products 
of the Z boson were taken from simulations with the |57|] 
Monte Carlo event generator using the tuned parameter 



set of the OPAL Collaboration |58 



The main ambiguities in propagation of photons con- 
cern the unknown rms magnetic field strength B which 
can influence the predicted 7— ray spectra via syn- 
chrotron cooling of the electrons in the EM cascade, and 
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the strength of the universal radio background (URB) 
which influences pair production by UHE 7— rays [j59| . 
Photon interactions in the GeV to TeV range are dom- 
inated by infrared and optical universal photon back- 
grounds (IR/O), for which we took the results of Ref. 
j60j. The resulting photon flux in GeV range in not sen- 
sitive to details of the IR/O backgrounds. 

Concerning the cosmology parameters we chose the 
Hubble parameter H = 70 km s _1 Mpc -1 and a cosmo- 
logical constant Qa = 0.7, as favored today. These values 
will be used in all cases unless otherwise indicated. 

For the neutrinos we assume for simplicity that all 
three flavors are completely mixed as suggested by exper- 
iments |6l| and thus have equal fluxes. For each flavor 
we sum fluxes of particles and antiparticles. 

Predictions for the all fluxes in both codes agree within 
tens of percents. Only for the Z-burst scenarios the pho- 
ton fluxes agree only within a factor ~ 2 between the two 
original codes due to the different implementation of Z- 
decay spectra and the ambiguities in photon propagation 
mentioned above. This difference has no influence on the 
conclusions of this paper. 

In the present investigation we parameterize power law 
injection spectra of either protons (for UHECR sources) 
or neutrinos (for Z-burst models) per co-moving volume 
in the following way: 

<P{E,z) = f(l + z) m E- a Q(E m ^~E) 

^min — Z — ^max ; (1) 

where / is the normalization that has to be fitted to 
the data. The free parameters are the spectral index 
a, the maximal energy £ mM , the minimal and maximal 
redshifts z m - m , z max , and the redshift evolution index to. 
The resulting neutrino spectra depend insignificantly on 
z nlin in the range < z min < 0.1 where local effects could 
play a role, and thus we will set z min = in the following. 

To obtain the maximal neutrino fluxes for a given set 
of values for all these parameters , we determine the max- 
imal normalization / in Eq. (|l|) by demanding that both 
the accompanying nucleon and 7— ray fluxes are below 
the observed cosmic ray spectrum and the diffuse 7— ray 
background observed by EGRET, respectively. 

III. THE COSMOGENIC NEUTRINO FLUX 
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FIG. 1: Dependence of the average maximal cosmogenic neu- 
trino flux per flavor consistent with all cosmic and 7— ray data 
on the maximal redshift z max , for the values indicated. Val- 
ues assumed for the other parameters in Eq. (jjj) for proton 
primaries are £ max = 10 23 eV, m = 3, a = 1.5. Also shown 
are the AGASA cosmic ray data above 3 x 10 18 eV Q. 
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FIG. 2: Dependence of the average maximal cosmogenic neu- 
trino flux per flavor consistent with all cosmic and 7— ray data 
on the maximal injection energy -E max , for the values indi- 
cated. Values assumed for the other parameters are z max = 2, 
a = 1.5, m — 3. The cosmic ray data are as in Fig. hi 



A. Dependence on unknown parameters 

In this section we discuss the case when primary UHE- 
CRs produce cosmogenic neutrinos as well as 7— rays 
during propagation. For EM propagation we use B = 
10 -9 G and the intermediate URB strength estimate of 
Ref. . These parameters only influence the 7— ray flux 
at UHEs, but not in the GeV range where the flux only 
depends on the total injected EM energy. Therefore, in 
this scenario the resulting neutrino fluxes are insensitive 
to the poorly known UHE 7— ray absorption because the 



"visible" UHE flux is always dominated by the primary 
cosmic rays and not by the secondary 7— ray flux, as can 
be seen in Figs. M and HO below. 

Figs. |], H ||, ancl^| show the dependencies of the 
cosmogenic neutrino flux (average per flavor) for which 
associated cosmic and 7— ray fluxes are consistent with 
the data, on the maximal source redshift z max , maxi- 
mal injection energy -E max , redshift evolution index to, 
spectral power law index a, and cosmological parame- 
ters, respectively. In each figure the line for the highest 
neutrino flux corresponds to a significant contribution of 
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FIG. 3: Dependence of the average maximal cosmogenic neu- 
trino flux per flavor consistent with all cosmic and 7— ray 
data on the source evolution index m, for the values indi- 
cated. Values assumed for the other parameters are z max = 2, 
E rnaJ r = 3 x 10 22 eV, a — 1. The cosmic ray data are as in 
Fig. @. 
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FIG. 4: Dependence of the average maximal cosmogenic neu- 
trino flux per flavor consistent with all cosmic and 7— ray 
data on the injection spectrum power law index a, for the 
values indicated. Values assumed for the other parameters 
are z max = 2, i5 max = 3 x 10 22 eV, m = 3. The cosmic ray 
data are as in Fig. |l[ 

the accompanying 7-ray flux to the observed flux at the 
EGRET region, whereas for the other lines the 7-ray flux 
gives negligible contributions to the EGRET flux. 

Fig. |2| shows that a change of the primary proton max- 
imum energy changes the secondary neutrino flux only 
in the high energy region. The reason for this behav- 
ior is the shape of the pion production cross section 
which at the lowest energies is dominated by the single 
pion A— resonance. Figs, [j] and || show that the cosmo- 
genic neutrino flux at the lowest energies mostly depends 
on the maximum redshift z max and the evolution index 
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FIG. 5: Dependence of the average maximal cosmogenic neu- 
trino flux per flavor consistent with all cosmic and 7— ray 
data on the cosmological vacuum energy density £7a and Hub- 
ble rate -ffo, measured in kms _1 Mpc _1 , for the values indi- 
cated. Values assumed for the other parameters are z max = 2, 
fimiix = 10 23 eV, q = 1.5, m = 3. The cosmic ray data are as 
in Fig. §. 

to. This is especially relevant for experiments with their 
main sensitivity below ~ 10 19 eV such as ICECUBE (see 
Fig. [h] below). Fig. |J shows that the maximal neutrino 
flux significantly increases with decreasing proton injec- 
tion power law index. Fig. [| shows that the variation 
of the maximal neutrino flux with cosmology parameters 
J7a and H is rather modest, about 50% for values dis- 
cussed in recent years. 



B. Active galactic nuclei as UHECR sources 

Here we consider AGN sources for the primary UHECR 
flux, with the typical evolution parameters m = 3.4 for 
z < 1.9 and to = for 1.9 < z < 2.9 §§. The only free 
remaining parameters are the power law index a and the 
maximum energy -E max for the proton injection spectrum, 
and the flux normalization / in Eq. (|l|) . We first checked 
that for the case a — 2 we agree with the WB bound. Our 
cosmogenic neutrino fluxes agree reasonably well with the 
corresponding ones shown in Fig. 4 of Ref. |33| when tak- 
ing into account maximal mixing assumed in the present 
paper. 

Fig. H shows the maximal cosmogenic neutrino flux for 
sources with hard spectra, a — 1, as a function of the 
maximal proton energy E max and a mono-energetic in- 
jection spectrum at E = 3 x 10 21 eV. For the case a = 1, 
MPR ||5q| computed the sum of the cosmogenic neutrino 
flux and the neutrino flux directly emitted by the sources 
which are assumed to be transparent to neutrons. Our 
fluxes shown here only include the cosmogenic flux. Nev- 
ertheless our fluxes overshoot the comparable MPR fluxes 
by up to a factor 5 at energies below the peak flux. This is 
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FIG. 6: Dependence of the average cosmogenic neutrino flux 
per flavor on maximum injection energy -E ma x, for the values 
indicated, assuming a — 1 and the AGN evolution parame- 
ters discussed in the text, "mono" indicates mono-energetic 
proton injection at E = 3 x 10 21 eV. For comparison, the 
7— ray bound derived from the EGRET GeV 7— ray flux , 
Eq.fcl), the MPR limit for optically thin sources |5q], and the 
WB limit for AGN-like redshift evolution S] are also shown. 
The cosmic ray data are as in Fig. hi. 



most likely due to a combination of the different cosmol- 
ogy (see Fig. |^) and a different implementation of multi- 
pion production which influences interactions of nucleons 
at high energies, thus at high redshift and in turn the low 
energy tail of cosmogenic neutrinos. 

Fig. ^ and also Figs. |^-[| demonstrate in a general way 
that it is easy to exceed the WB bound and even the 
MPR bound for injection spectra harder than about E~ 2 . 
This is because Waxman-Bahcall restricted themselves 
to nucleon injection spectra softer than E~ 2 and sources 
smaller than nucleon interaction lengths |55| . Thus, their 
bound does not directly apply to the cosmogenic neu- 
trino flux. In addition, in our opinion, their assumptions 
on the injection spectra are too narrow: Possible sce- 
narios with hard injection spectra and the AGN redshift 
evolution assumed here (which is the same as the one 
used by Waxman-Bahcall) include cases where UHECRs 
are accelerated by the electromotive force produced by 
magnetic fields threading the horizons of spinning super- 
massive black holes in the centers of galaxies [Q or by 
reconnection events around forming galaxies ]65|. 

In any scenario involving pion production for the cre- 
ation of 7— rays and neutrinos, the fluxes per flavor are 
approximately related by F V {E) ~ F 1 (E)/3. Assuming 
smooth spectra and comparing with the EGRET 7— ray 
fluence, energy conservation implies 



E 2 F V (E) < 6 x 10 2 eVcm 



(2) 



This ultimate bound is also shown in Fig. ^|and Fig. |7 be- 
low. It corresponds to the MPR limit for optically thick 
The maximal E 2 j(E) of the fluxes in Figs. @ 



and [?] below indeed reach this 7— ray bound Eq. (||). 

Note that the two theoretical bounds shown in Fig. || 
represent fluxes per neutrino flavor under the assump- 
tion of maximal neutrino mixing. They are thus about 
a factor 2 lower than in Refs. |^5|, |56) which show 
muon neutrino fluxes in the absence of mixing where the 
tau-neutrinos fluxes are negligible and electron neutrino 
fluxes are about a factor 2 smaller. The WB and MPR 
bounds represent upper neutrino flux limits for compact 
sources such as AGN and 7— ray bursts in case of small 
optical depth for nucleons. We discuss a specific non- 
shock acceleration scenario in Sect. VII, for which both 
of the above bounds are not valid because the source 
optical depth for nucleons is large (for reviews on AGN 
neutrino fluxes see, e.g., Ref. (|]). 



C. General case of arbitrary source evolution 
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FIG. 7: Dependence of the average cosmogenic neutrino flux 
per flavor maximized over maximal injection energy -E max , 
evolution index m, and normalization consistent with all cos- 
mic and 7— ray data, on the injection spectrum power law 
index a. "mono" indicates mono-energetic proton injection 
at E = 10 21 eV. The rest of the line key is as in Fig. [| 

In this section we consider more general UHECR 
sources and relax the restrictions on their redshift evolu- 
tion. Fig. ^ shows that cosmogenic neutrino fluxes higher 
than both the WB and MPR limits are possible even for 
relatively soft E~ 2 proton injection spectra, if the red- 
shift evolution is stronger than for AGNs: The curve 
for E~ 2 in Fig. |^ corresponds to the evolution parame- 
ters m = 5, z mm = 3 and E max = 10 22 eV, the curve 
for E- 1 - 75 to rn = 4.5, z max = 3 and E max = 10 23 eV, 
and the curve for mono-energetic injection to m = 4, 
W = 3, and E max = 10 21 eV. 

Fig. I shows that between ~ 10 18 eV and ~ 10 20 eV the 
energy loss rate of protons on the CMB is dominated by 
pair production instead of pion production. The former 
does not contribute to neutrino production but the EM 
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FIG. 8: The nucleon interaction length (dashed line) and en- 
ergy attenuation length (solid line) for photo-pion production 
and the proton attenuation length for pair production (thin 
solid line) in the combined CMB and the estimated total ex- 
tragalactic radio background intensity. 



cascades initiated by the pairs lead to contributions to 
the diffuse 7— ray background in the GeV range. Thus, 
the cosmogenic neutrino flux is the more severely con- 
strained the bigger the fraction of cosmic ray power is 
in the range 10 18 eV < E < 10 20 eV. This is mostly im- 
portant for soft injection spectra and explains why the 
total neutrino energy fluence decreases with increasing a 
in Fig. @. 



D. Comparison with experimental limits and 
future sensitivities 



Figs. [5] and [l^ shows a scenario maximized over all 4 
parameters in comparison to existing neutrino flux up- 
per limits and expected sensitivities of future projects, 
respectively. Both of these fall into two groups, detec- 
tion in water or ice or underground, typically sensitive 
below ~ 10 16 eV, and air shower detection, usually sen- 
sitive at higher energies. Existing upper limits come 
from the underground MACRO experiment jllj at Gran 
basso, AMANDA [| in the South Pole ice, and the Lake 
BAIKAL neutrino telescope in the first category, and 
the AGASA ground array j30| , the former fluorescence 
experiment Fly's Eye J2{|, the Radio Ice Cerenkov Ex- 
periment RICE j}4|, and the Goldstonc Lunar Ultra- 
high energy neutrino experiment GLUE [ {35| in the sec- 
ond category. Future experiments in the first category 
include NT200+ at Lake Baikal §, ANTARES in the 
Mediterranean [§, NESTOR in Greece @, AMANDA- 
II [[70), and ICECUBE ||, the next-generation ver- 
sion of AMANDA, whereas the air shower based cat- 
egory includes the Auger project p8|, the Japanese 
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7), the fluorescence/Cerenkov detec- 
and the space based OWL |38| and 



FIG. 9: A scenario with maximal cosmogenic neutrino fluxes 
as obtained by tuning the parameters to z max = 2, i5 max = 
10 23 eV, m = 3, a = 1. Also shown are predicted and ob- 
served cosmic ray and 7— ray fluxes, the atmospheric neutrino 
flux j(5(|, as well as existing up per limits on the diffuse neu- 
trino fluxes from MACRO jul AMANDA 0L BAIKAL §, 
AGASA @, the Fly's Eye |S) and RICE [MJ experiments, 
and the limit obtained with the Goldstone radio telescope 
(GLUE) Pq |, as indicated. The cosmic ray data are as in 
Fig. I 



EUSO |24j experiments. The vertical bars for the 
MOUNT sensitivity characterize the uncertainties due 
to the not yet determined zenith angle range and sen- 
sitivity to the fluorescence component. The OWL sen- 
sitivity estimate is based on deeply penetrating atmo- 
spheric showers induced by electron or muon-neutrinos 
only |58| and may thus be considerably better if tau 
neutrinos, Cerenkov events, and Earth skimming events 
are taken into account J69| . The same applies to the 
EUSO project ||. The AMANDA-II sensitivity will lie 
somewhere between the ANTARES and ICECUBE sen- 
sitivities [fzo|| - The maximized fluxes shown in Figs. ^ 
and [h] arc considerably higher than the ones discussed 
in Refs. |37], fn], [72| [73), and should be easily de- 
tectable by at least some of these future instruments, as 
is obvious from Fig. RG. 



IV. NEUTRINO FLUXES IN TOP-DOWN 
SCENARIOS 

Historically, top-down scenarios were proposed as an 
alternative to acceleration scenarios to explain the huge 
energies up to 3x 10 20 eV observed in the cosmic ray spec- 
trum j74| . In these top-down scenarios UHECRs are the 
decay products of some super-massive "X" particles of 
mass mx ^ 10 20 eV close to the grand unified scale, and 
have energies all the way up to ~ mx- Thus, The mas- 
sive X particles could be metastable relics of the early 
Universe with lifetimes of the order the current age of 
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FIG. 10: The cosmogenic neutrino flux (solid line) shown in 
Fig. ^ in comparison with expected sensitivities ol the cur- 
rently being constructed Auger project to electron/muon and 
tau- neutrinos pq |, and the planned projects telescope array 
(TA) |j7| (dashed-dotted line), the fluorescence/Cerenkov de- 
tector MOUNT [32) . and, indicated by squares, the space 
based OWL Jij]] (we take the latter as representative also 
for EUSO).the water-based NT200+ §, ANTARES [§ (the 
NESTOR @ sensitivity would be similar to ANTARES ac- 
cording to Rel. 0), and the ice-based ICECUBE ||, as indi- 
cated. Also shown (dashed line) is an extreme scenario with 
Zmax = 3, Umax = 10 22 eV, 77i = 5, and a — 2, leading to a 
cosmogenic neutrino flux extending to relatively low energies 
where ANTARES and ICECUBE will be sensitive, and the 
atmospheric neutrino flux for comparison. 



the Universe or could be released from topological de- 
fects that were produced in the early Universe during 
symmetry-breaking phase transitions envisaged in Grand 
Unified Theories (GUTs). The X particles typically de- 
cay into leptons and quarks. The quarks hadronize, pro- 
ducing jets of hadrons which, together with the decay 
products of the unstable leptons, result in a large cas- 
cade of energetic photons, neutrinos and light leptons 
with a small fraction of protons and neutrons, some of 
which contribute to the observed UHECR flux. The re- 
sulting injection spectra have been calculated from QCD 
in various approximations, see Ref. JT^ ] for a review and 
Rcf. ffdji for more recent work. In the present work we 
will use the spectra shown in Fig. [ll| which are obtained 
from solving the DGLAP equations in modified leading 
logarithmic approximation (MLLA) without supersym- 
metry for X particles decaying into two quarks, assuming 
10% nucleons in the fragmentation spectrum. 

For dimensional reasons the spatially averaged X par- 
ticle injection rate can only depend on the mass scale mx 
and on cosmic time t in the combination 



hx(t) = Km p x t 



-i+P 



(3) 



FIG. 11: Unnormalized nucleon, electron, 7— ray, and neu- 
trino (per flavor) MLLA spectra resulting from X particles de- 
caying into two quarks without supersymmetry. These spec- 
tra are used as injection spectra of top-down models in the 
present work. The spectra of antiparticles can be assumed 
identical. 



whereas decaying super-heavy dark matter of lifetime 
much larger than the age of the Universe corresponds 
to p = 2 |Q . In the latter case the observable flux will 
be dominated by decaying particles in the galactic halo 
and thus at distances smaller than all relevant interaction 
lengths. Composition and spectra will thus be directly 
given by the injection spectra which are most likely in- 
consistent with upper limits on the UHE photon fraction 
above 10 19 eV |7(| , see Fig. |nj. In addition, decaying dark 
matter scenarios suffer in general from a more severe fine 
tuning problem and predict a smaller neutrino flux than 
extragalactic topological defect model scenarios. We will 
therefore here focus on the latter, with p = 1. 

Fig. |l^ shows the results for m x — 2 x 10 14 GeV, with 
B = 1CU 12 G, and again the minimal URB consistent 
with data |59|, [77|. These parameters lead to optimistic 
neutrino fluxes for the maximal normalization consistent 
with all data. For more detailed recent discussions of 
top-down fluxes see Refs. jr|], (79|. 



V. THE Z-BURST SCENARIO WITH 
ACCELERATION SOURCES 

In the Z-burst scenario UHECRs are produced by 
Z-bosons decaying within the distance relevant for the 
GZK effect. These Z-bosons are in turn produced by 
UHE neutrinos interacting with the relic neutrino back- 
If the relic neutrinos have a mass m„, Z- 



ground |45 



where k and p are dimensionless constants whose value 
depend on the specific top-down scenario |74|]. Extra- 
galactic topological defect sources usually predict p = 1, 



bosons can be resonantly produced by UHE neutrinos 
of energy E v ~ M%/{2m v ) ~ 4.2 x 10 21 eV (eV/m„). 
The required neutrino beams could be produced as sec- 
ondaries of protons accelerated in high-redshift sources. 
The fluxes predicted in these scenarios have recently been 
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FIG. 12: Flux predictions for a TD model characterized by 
P = j-' m x = 2 x 10 14 GeV, and the injection spectra given in 
Fig. [H]. The line key is as in Fig. 



discussed in detail in Refs. []46|, |8CJ] . In Fig. |T3| we show 
an optimistic example taken from Rcf. po| for compar- 
ison with the other scenarios. As in Refs. jl6| |8(| no 
local neutrino over-density was assumed. The sources 
are are assumed to not emit any 7— rays, otherwise the 
Z-burst model with acceleration scenarios is ruled out, as 

We note that no known astro- 



demonstrated in Ref . 1 46 



physical accelerator exists that meets the requirements of 
the Z-burst model G3]. Also note that even exclusively 
emitting neutrino sources appear close to be ruled out 
already by the GLUE bound. 



FIG. 13: Flux predictions for a Z-burst model averaged over 
flavors and characterized by the injection parameters z m i n = 
0, Zmax = 3, a = 1, m = 0, in Eq. ([!]) for neutrino primaries. 
All neutrino masses were assumed equal with m v = 0.1 eV 
and we again assumed maximal mixing between all flavors. 
The line key is as in Fig. bl 



to a decrease of the UHE proton and photon fluxes, see 
Fig. [l4|. The EGRET measurement therefore consider- 
ably constrains the parameter space of this model. Neu- 
trino masses m„ 3> 0.1 eV are required, which allow 
X-particle masses mx % 10 14 GeV implying secondary 
photon fluxes that are below the measured level in the 
EGRET energy range. 



VI. MONO-ENERGETIC SUPER-HEAVY RELIC 
NEUTRINO SOURCES 

Since no known astrophysical accelerator exists that 
produces sufficiently strong neutrino beams up to suffi- 
ciently high energies for the Z-burst scenario to work [^6[ , 
one can speculate about more exotic sources. Gelmini 
and Kusenko jl7| have considered a top-down type source 
for the Z-burst scenario, namely super-heavy particles 
mono-energetically decaying exclusively into neutrinos. 
In Fig. [l4| we show predictions for one example of this 
kind of model with the same maximal energy -E max = 
10 23 eV = rnx/2 as for case of Z-burst model. Again, 
no local neutrino over-density was assumed and the cal- 
culation assumed a minimal URB flux and a magnetic 
field B = 10~ 12 G. Note that, for the same UHECR flux, 
in this scenario the photon flux in the EGRET region is 
larger than in the Z-burst model with power-law neutrino 
sources, compare Fig. [l3[ because all secondary photons 
from Z-boson decays at redshifts z > 3 contribute only 
to the EGRET energy range. Thus, the normalization 
of the photon flux to the EGRET measurements leads 




E[eV] 

FIG. 14: Flux predictions for a Gelmini-Kusenko model char- 
acterized by p = 2, mx = 2 x 10 14 GeV in Eq. (|), with 
X particles exclusively decaying into neutrino-anti-neutrino 
pairs of all flavors (with equal branching ratio), assuming all 
neutrino masses m v =0.1 eV. The line key is as in Fig. |9[ 

Let us also note though the neutrino flux in UHECR 
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region is much smaller in Gclmini-Kusenko model in com- 
pare with Z-burst model, the GLUE bound constrains 
both models in the peak of neutrino flux in similar way. 



VII. NEUTRINO FLUXES IN A NON-SHOCK 
ACCELERATION MODEL FOR AGN 
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FIG. 15: Neutrino flux predictions for the AGN model H for 
a uniform distribution of blazars (no redshift evolution )~The 
photon flux is below measured EGRET value. The typical 
neutrino flux in this model contains the same energy as the 
photons. The position of the peak is governed by the initial 
proton distribution. The line key is as in Fig. [To| 

Although the recent exciting discoveries by the Chan- 
dra X-ray space observatory added much to our knowl- 
edge of structures of the jets of powerful radiogalaxies 
and quasars, they didn't solve the old problems and, in 
fact, brought new puzzles. If the observed X-ray emission 
is due to synchrotron energy losses of electrons, the en- 
ergy of such electrons should be of the order of 100 TeV. 
Electrons of such energies loose all their energy on a typ- 
ical scale of only 0.1 kpc. In order to explain observed 
X-ray data, such 100 TeV electrons should be created 
uniformly over the jet length of order of 100 kpc. The 
conventional shock-wave acceleration mechanism in this 
case would require a jet uniformly filled with 1000 iden- 
tical shocks, following one another along the jet axis. 

In Rcf. |Q a "non-shock acceleration" version of the 
electron synchrotron model was proposed, namely assum- 
ing that electrons are not accelerated in the jet, but are 
instead the result of pair production by very high energy 
(VHE) 7— rays interacting with the CMB. The typical 
attenuation length of 10 16-18 eV photons is of order 100 
kpc, comparable with the lengths of large scale extra- 
galactic jets. 

In this model the VHE 7-rays are produced by acceler- 
ated protons interacting with the ambient photon fields 



(supplied, for example by the accretion disk around the 
massive black hole) through photo-meson processes. At 
the same time those protons produce neutrinos which are 
emitted in the direction of the jet. Therefore, this model 
predicts a high neutrino flux comparable in power with 
the 7-ray flux. The detailed numerical simulations of 
proton acceleration in the central engine of the AGN 
show that the collimated jet of almost mono-energetic 
VHE protons (linear accelerator) can be created in the 
electro-magnetic field around the black hole and the en- 
ergy of those protons can be converted into photons and 
neutrinos, while protons can be captured inside of the 
source. The nucleon flux leaving the AGN is well below 
observed cosmic ray flux in this scenario. Furthermore, 
since all nucleons leaving the source are well below the 
GZK cutoff, there is no cosmogenic contribution to the 
neutrino flux from these sources. 

Fig. |l^ shows a typical prediction of the diffuse neu- 
trino flux in this scenario. This flux is beyond the WB 
limit which is not applicable in this case because the 
sources are optically thick for nucleons with respect to 
pion photo-production. The flux is consistent with MPR 
bound for optically thick sources. 

In the AGN model discussed above, blazars would be 
seen by neutrino telescopes as point-like sources with 
neutrino fluxes which are smaller or of the same order 
as the photon flux emitted by these same sources and 
which are detectable by 7-ray telescopes. 



VIII. CONCLUSIONS 

Based on our transport code we reconsidered neutrino 
flux predictions and especially their maxima consistent 
with all cosmic and 7— ray data, for cosmogenic neutrinos 
produced through pion production of UHECRs during 
propagation, and for the more speculative Z-burst sce- 
nario and top-down scenarios. We pointed out that one 
can easily exceed the WB bound and, in the most opti- 
mistic cases, even the MPR bound for cosmogenic neutri- 
nos in scenarios with cosmic ray injection spectra harder 
than E^ 1 - 5 , maximal energies fi max > 10 22 eV, and red- 
shift evolution typical for quasars, or stronger. Given our 
poor knowledge on the origin of UHECRs, in our opin- 
ion these are possibilities that should not be discarded at 
present, especially since they would lead to considerably 
increased prospects of ultra-high energy neutrino detec- 
tion in the near future. We also show that for non-shock 
AGN acceleration models the AGN neutrino fluxes can 
reach the 7— ray bound Eq. (0) around 10 16 eV which 
represents the ultimate limit for all scenarios of 7— ray 
and neutrino production involving pion production. 

Finally, we note that fluxes as high as for the optimistic 
scenarios discussed here would also lead to much stronger 
constraints on the neutrino-nucleon cross section at ener- 
gies beyond the electroweak scale. This is important, for 
example, in the context of theories with extra dimensions 
and a fundamental gravity scale in the TeV range p^|. 
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